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ABSTRACT 

A radial alignment study of three solar wind stream structures observed by 
IMP-7 and -8 (at t.O AU) and Voyager 1 and 2 (in the range 1.4 - 1,8 AU) in 
late 1977 is presented. It is demonstrated that several important aspects of 
the observed dynamical evolution can be both qualitatively and quantitatively 
described with a single-fluid 2-D MHD numerical model of quasi-steady 
corotating flow, including accurate prediction of: 1) the formation of a 

corotating shock pair at 1.75 AU in the case of a simple, quasi-steady stream; 
2) the coalescence of the thermodynamic and magnetic structures associated 
with the compression regions of two neighboring, interacting, corotating 
streams; and 3) the dynamical "destruction" of a "small" (i.e., low 
velocity-amplitude, short spatial-scale) stream by its overtaking of a slower 
moving, high-density region associated with a preceding transient flow. The 
evolution of these flow systems is discussed in terms of the concepts of 
"filtering" and "entrainment", which have appeared in the literature. It is 
concluded that this study provides support at a detailed, quantitative level 
for our general conceptual notions of stream evolution in the inner solar 
system. 
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1. Introduction 


The 3tudy of dynamics of interplanetary stream evolution has progressed to 
the point that it can by legitimately claimed that the basic processes are 
fairly well understood. It was early appreciated (Parker, 1963; Sarabhai, 
1963) that large scale inhomogeneities in solar wind flow speed associated 
with the quasi-steady global structure of the corona would, under the 
influence of solar rotation, give rise to a spiral interaction front in 
interplanetary space. At the leading edge of the structure, kinematic 
steepening would occur as fast material overtakes slower material previously 
emitted from the corona, thereby compressing the gas caught between the two 
flows. Likewise, at the trailing edge of the high-speed scream, a rarefaction 
would form as fast material outruns the succeeding slow flow. If the speed 
differences in the two flows were great enough, forward and reverse shock 
waves could be expected to form at the eri^es of the compression region. 
Quantitative models (e.g., Hundhausen, 1973a) suggested shock pairs might form 
as near to the sun as 1 AU, and simple dynamical considerations dictated these 
structures should survive well out into the solar system (Hundhausen, 1973b). 

A reverse shock was observed at 1 AU by Burlaga, 1979, but in general, 
corotating shocks are rarely seen at 1 AU (Ogilvie, 1972). Observations by 
Pioneer 10 and 11 subsequently confirmed that shock pairs form between 1 AU 
and 4 or 5 AU (Hundhausen and Gosling, 1976; Smith and Wolfe, 1976, 1977; 

Dryer et al., 1978), where numerous sawtooth velocity structures were 
detected. These structures (dubbed corotating interaction regions or CIR's) 
represent the more advanced stages of stream evolution, in which the material 
entrapped at the stream front has achieved its maximum compression and is 
reexpanding, driving a forward-reverse shock pair back into the surrounding 
flow. It was found the entire process could be nicely modeled with a 
relatively simple nonlinear hydrodynamic formulation (Gosling et al., 1976), 
though some approximation to the field behavior would be obviously required to 
provide a more complete analysis (e.g., Dryer et al., 1978). 

Though the gross aspects of stream evolution appear to have been 
satisfactorily explained, a great many details remain to be worked out. In 
recent years, for example, considerable attention has been devoted to the 
quantitative evaluation of the dynamic mechanisms by which the gas at the 
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leading edge of a high-speed stream responds to the kinematic steepening 
imposed by the velocity structure (Goldstei", 1971* Pizzo, 1980, 1981, 1982; 
Whang and Chien, 1981). Gas and magnetic pressure gradients built up in the 
compression region give rise to nonradial flows, which leads to lateral 
slippage between the colliding fast and slow streams. The net effect i3 that 
the dynamic reaction of the gas tends to oppose the kinematic steepening, 
resulting in a less severe, spatially broadened compression and a 
correspondingly reduced rate of shock formation. These effects are thought to 
be particularly important out to about 1 AU and at least partially explain 
some of the not insubstantial discrepancies between theory and observation in 
that region. (See Pizzo, 1983, for a more complete discussion of this 
problem.) 

In view of the wealth of data and more sophisticated models currently 
available, it would be instructive to combine the two in taking a harder and 
more detailed look at the various dynamical processes occurring in the 
interplanetary solar wind. To this end, we have assembled a set of plasma and 
field data from 1977, when Voyager 1 and 2 were situated between 1.4 and 1.8 
AU and were favorably radially aligned with the earth-orbiting IMP-7 and -8 
spacecraft as well as with each other. Three successive stream structures 
from this period will be analyzed with the help of a single fluid, 2-D MHD 
corotating flow model to illustrate several interesting features of the 
evolution and to provide a theoretical basis for their interpretation. 

The first structure to be discussed is the simplest one, appearing at 1 AU 
as a classic steepened (but not shocked) quasi-steady high speed stream. It 
is demonstrated that the model successfully predicts the formation of a 
forward-reverse shock pair that was subsequently observed by Voyager 1 at 1.74 
AU and also reproduces the observed plasma and magnetic variations with 
reasonable accuracy. 

The second case involves a more complicated situation, in which one fast 
corotating stream overtakes a slower, preceding corotating stream. While the 
same physics operative in the simple stream still applies, it should come as 
no surprise that certain details of the process will differ. Here both the 
faster and slower streams are individually evolving at somewhat different 
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rates* while at the same time the faster, ’’bigger" one advances upon (and 
sweeps up) the slower, "smaller" one. We shall be paying particular attention 
to the thermodynamic and magnetic structures of the two streams, for as the 
evolution proceeds the compression regions of the individual streams are 
eventually advected together, resulting in a single, merged compression at the 
front of a single, merged stream. Again, the model is shown to describe the 
observed processes quite well. 

In the third case, we consider what happens when a "small" corotating 
stream overtakes a transient high density region behind a shock. Comparison 
of the observations with the model demonstrates that the stream is effectively 
destroyed in the interaction, i.e,, that all the high speed material in the 
stream is consumed by the dynamic exchange of momentum with the dense material 
of the transient flow. 

Together, these studies bolster our confidence in the ability of 
2-dimensional MHD models to replicate the dynamic behavior of the solar wind 
near and beyond 1 AU and lend further credence to our general conceptual 
picture of stream evolution. Moreover, they also serve as excellent 
illustrations of the dynamical processes that have been variously described in 
the literature as either "filtering" (Gosling et al., 1976) or "entrainment" 
(Burlage et al., 1983a) and will therefore be used as a springboard tv a 
discussion of this subject in the final passages of this paper. 



2 • Overview of the Flows, Spacecraft Trajectories and Data Sets 


For the purposes of this study, we shall make use of the solar wind flows 
observed by IMP 7 and 8 at 1 AU and Voyager 1 and 2 at 1. <1-1.7 AU in the time 
period between November 10 and * December 25, 1977. The speed profile 
observed by IMP-7 and -8 at 1 AU in this period is plotted in Figure 1 
together with a map of coronal holes from the HAO. The coronal hole map has 
been shifted in time corresponding to the transit time to 1 AU, so that the 
streams and the corresponding coronal holes are nearly vertically aligned. 
Three flows predominate in the interval. The first and third flows (labeled 1 
and 3 in Figure 1) apparently originated in an equatorial extension of a lobe 
of the north polar coronal hole, and they represent a recurrent and hence a 
"corotating flow". Note, however, that neither the streams nor the coronal 
holes are unchanged from one rotation to the next), through they are steady 
enough to be classified as corotating flows, on the basis of solar wind data. 
Such flows were discussed by Burlaga et al. (1978) and Rhodes and Smith 
(1981). Stream 2 is associated with a lobe of the south polar coronal hole, 
and it has negative polarity for the most part, 

Closer inspection of the magnetic field and speed profiles of streams 1 
and 2 suggests that they are actually "compound streams" (Burlaga and Ogilvie, 
1970), each consisting of two separate streams. The initial part of stream 1 
has "anomalous" negative polarity, whereas the body of the stream has positive 
polarity consistent with its source being a positive polarity, northern 
hemisphere coronal hole. Similarly, the initial part of stream 2 has mixed 
polarity, whereas the region with highest speeds has negative polarity 
consistent with its source being a southern hemisphere coronal hole. 
(Cross-hatched regions in Figure 1 indicate missing magnetic field data.) 
Stream 3 is relatively uncomplicated, with uniform positive polarity, and it 
can be thought of as a simple corotating stream. 

The trajectories of Voyagers 1 and 2 (hereafter abbreviated VI and V2, 
respectively) and of IMPs 7 and 8 are shown in inertial heliographic (IHG) 
coordinates in Figure 2. X TM _ and Y_„_ are the coordinate axes in the solar 
equatorial plane: X points along the longitude of the ascending node. The 

latitude angle \ between the Sun-spacecraft line and the solar equatorial 
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plane is plotted versus X IHG in Figure 2. Owing to the tilt of the sun’s 
equatorial plane relative to the ecliptic, the heliographic latitude of Earth 
(hence also IMP-7 and IMP-8) oscillates between +7.5° and -7.5° in the course 
of the year, VI was very close in latitude to IMP-7 and IMP-8 throughout most 
of the interval of interest, while V2 was at somewhat higher latitudes. IMP-7 
and IMP-8 were radially aligned with VI and V2 near November 20, and the 
alignment was close to radial throughout the time interval, the misalignment 
being greatest for stream 3, Since the IMP-V1 alignment was generally better 
and since there are only minor differences between the VI and V2 observations, 
for clarity only the VI data will be used in the plots. 

This study was based primarily upon hourly averages of the magnetic field 
strength and components from the GSFC magnetometers on IMP-8, VI and V2, and 
hourly averages of the proton density, temperature and bulk velocity from the 
MIT plasma experiments on IMP-7, IMP-8, VI and V2; some higher resolution data 
will be discussed in relation to the identificataion of shocks and 
discontinuities. Owing to instrumental problems, magnetic field data and 
measurements of density and temperature are not available from IMP-7, so there 
are gaps of *4 days in those parameters every ^8 days (when IMP-8 was in the 
magnetosphere). Thus, our 1 AU baseline data consist of nearly continuous 
measurements of the bulk speed from IMPs-7 and -8, and measurements of the 
magnetic field, density, and proton temperature from IMP-8 alone. There are 
also small gaps of less than a few hours duration every day or so, but these 
can reasonably be filled by linear extrapolation, since the stream profiles 
generally change very little on a scale of a few hours. These data are from 
the compilation of King (1979). The Voyager data set is more nearly complete, 
with measurements available nearly every day, but there are intermittent 
tracking gaps of several hours duration. Again these can generally be filled 
by extrapolation. 


3. The Model and Its Application to the Data 

An essential part of the analysis is to feed data from one spacecraft at 
lesser heliocentric distance into a dynamical model and map the evolution to a 
suitably aligned spacecraft at a greater heliocentric distance, and to compare 
the theoretical predictions with the second set of observations. Burlaga and 
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Ogilvie (1970) showed that the magnetic pressure is comparable to the thermal 
pressure in interaction regions at 1 AU, and they stressed that a 
magnetohydrodynamic theory is needed to describe their behavior. MHD models 
of flows beyond 1 AD have been presented by Dryer et al. (1976* 1978)* 
Goldstein and Jokipii (1977)* and Pizzo (1983). The numerical model employed 
here is a planar version of the single-fluid 3-0 formulation of Pizzo (1982). 
Since the model is described in great detail in the work oited, only a brief 
outline is presented here. 

The 2-D nonlinear, polytropic MHD equations are used to describe 
large-scale* quasi-steady, corotating flow dynamics in interplanetary space. 
In an inertial (non-rotating) frame, the single-fluid equations in 
conservation form are: 
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(Maxwell's equation) 
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where v is the azimuthal velocity in the rotating frame. 

The two independent variables are the spherical components (r,$), and the 

five dependent variables are the radial and azimuthal velocity components u r 

and u., the radial magnetic field B , and the two thermodynamic variables, 

<P r 

density p and gas pressure P. Pressure is defined by the perfect gas law as 

P s 2nkT, (6) 

where n is the proton number density (= p/m ), k is Boltzmann's constant, and 

P 

T is the single- fluid temperature. For purposes of computation, T at 1 AU is 

obtained by taking the average of the observed proton temperature T and a 

4o p 

presumed constant electron temperature T =1.25 x 10 K. 
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Since the field is effectively tied to the flow in the rotating f^ame due 
the high electrical conductivity of the medium, the azimuthal component of the 
field is given by 


B = — ^ B (7) 

» u r 

In the 2-D approximation, only flow within the equatorial plane is considered, 
50 U Q s s 3,/30 = °* wh ' sre 9 is the meridional (north-south) coordinate in 
the spherical system. Thus all structures are assumed to be uniform along 
meridional arcs and are oriented perpendicular to the equatorial plane. 
Finally, the internal energy is dictated by a polytrope law [incorporated in 
the total energy equation (5)3. The polytropic index, y» for the single fluid 
gas is assigned the value 5/3; the results are not sensitive to the choice of 
Y, Two fluid models were alet: examined, but there was little difference 
between the predictions of this model and those of the 1-fluid model (see the 
Summary and Discussion). The remaining physical constants in ( 1 )— (5 ) are fl, 
the mean solar rotation rate; G, the universal gravitational constant; and M, 
the solar mass. 

The numerical Integration of the equations is fully discussed in Pizzo 
(1982) and only a brief synopsis is needed here. Values for the dependent 
variables are specified at each discrete $ at some initial r and are 
propagated to r + ar by means of the two-step explicit scheme of MacCormack 
(1971). The radial step size ar is determined from the Courant condition, 
taking into account additional constraints imposed by the small amounts of 
artificial viscosity introduced to maintain computational stability in the 
presence of corotating shocks and discontinuities. 

The incorporation of real data into the model requires a string of input 
data evenly spaced in azimuth (longitude), all gathered at a single, fixed 
heliocentric distance r. Hourly averages from the Earth-orbiting IMP 
spacecraft fit this prescription quite closely, requiring only a simple 
conversion from time to longitude. For the large-scale structures to be 
modeled here, the resulting resolution of about 0.5° is quite adequate to 
follow the evolution. Hourly averages from th** more distant Voyager 
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spacecraft provide the same azimuthal resolution, but an added complication 
arises from the radial motion of the spacecraft. The Voyagers typically 
travel one-tenth AU in heliocentric distance during the time it takes for the 
complete stream structure to pass by. For exactitude in the model-data 
comparison, we should either interpolate the model output to the instantaneous 
position of Voyager or, less palatably, correot the data to some common r by 
arbitary means. For the level of accuracy required in this study, though, we 
have elected to simply overlay on the raw Voyager time sequence the output of 
the model at a fixed r chosen such that the Voyager and model radii ooincide 
in the middle of the interaction region. In thi3 way, the largest radial 
position displacement between the two curves is limited to about .01 AU in the 
vicinity of the dynamic structures of greatest interest, The point where the 
data and model coordinates exactly coincide is indicated in the plots by the 
triangle along the bottom. 

Data gaps in the input have been filled by linear interpolation. Where 
only a few hours data are missing, little degradation in the mapping is 
experienced. However, where the gaps are more serious, particularly in the 
IMP record, spurious behavior is evident in the theoretical predictions. 
Fortunately, since the flow is highly supersonic and the radial span of the 
mapping is fairly small, these artifacts are relatively localized and the 
tainted portions of the solution can be readily identified and discarded. 
Similarly, the treatment of the longitudinal ends of the data string, which 
are arbitrarily forced to be periodic for computational considerations, has 
has negligible influence on the main part of the solution. 

In specifying the magnetic field input, we are presented with a choice, 
since the spacecraft data contain more information that is compatible with the 
corotating formulation. That is, once u and u. are set, then the 

P $ 

field-aligned condition (7) dictates that only one element of information 

about the field can be entered into the initial conditions. We could specify 

B or B directly from the data, but the behavior of the components of § in 
r 9 

the solar wind is generally dominated by waves and small-scale transient 
structures. This problem becomes especially acute at large heliocentric 
distances, where hourly averages of the components are not representative of 
the mean field and determination of the field magnitude at the lowest levels 
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might bo subjeot to uncertainty. On the other hand, it is |$| which is the 

important physical parameter, for the main large-scale magnetio driving foroe 

in interplanetary space is the magnetic pressure force, not tension. 

Furthermore, |Sj regulates the characteristic speed at which small-amplitude 

disturbances propagate and are a prime factor in the relaxation of CIR's* In 

the model calculations, on the initial surfaoe is fixed by the relation 

r 


B 


r 



, 2 
1 + -*5 

u r 


1/2 


r 


where |8|, v., and u are the observed quantities. Since all magnetic force 
9 r 

terms in the corotating formulation are independent of the polarity of B^, the 
positive sign is always chosen. Magnetic sectors may be passively mapped by 
following the motion of select tracer particles in the flow and reversing the 
polarity across the appropriate boundaries after the faot. Of course, the 
dyns-'iins associated with the current sheet transition are not accounted for, 
bo/* t vi use effects should be quite localized and negligible on the large scale. 
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4. Classic Corotating Stream 

We first consider the radial evolution of the simple corotating stream 
labeled "stream 3" in Figure 1, Five-minute averages of the basic parameters 
measured st 1 AU (Figure 3) show a stream interface and an interaction region 
ahead of a fast stream, but there is no indication of a forward or reverse 
shock. Hourly averages of these data, shown on the left of Figure 4, were 
used as inputs to the MHD model, and the predicted profiles at 2 AU are shown 
at the right of Figure 4. Note that the interface survives to 2 AU, and more 
importantly a forward and reverse shock pair are present at 2 AU, forming the 
boundaries of a relatively simple pressure wave. Voyager 1 observed this 
stream at 1.74 AU, and the observations as well as the MHD model predictions 
at this distance are shown in Figure 5. This figure shows that a corotating 
forward and reverse shock pair are both predicted and observed. The predicted 
shocks arrive slightly earlier than the observed shocks (an error of < 5° in 
Carrington longitude or < 10}( error in the propagation time between 1 AU and 
1,74 AU). Tills small difference can be attributed to variations associated 
with the departures from radial alignment and possibly temporal and 
latitudinal variations. 

The existence of the forward-reverse shock pair at 1.74 AU is shown more 
clearly in Figure 6 which gives a plot of 9.6s averages of the magnetic field 
and measurements of the plasma parameters. The forward shock is very thin, 
the transition occurring in < 9.6 sec. The reverse shock is somewhat broader, 
with the density and thermal speed decreasing over an interval of several 
minutes, but it is still appropriate to regard it as a shock. 

Returning to Figure 5, it is useful to compare the predicted and observed 
variations of some other flow parameters, the density, field strength, 
temperature, and total pressure. The longitudinal shift of the major features 
near the stream front is readily apparent, but the shapes of the profiles and 
their approximate amplitudes match reasonably well. A stream interface 
(marked I in the plot) is both predicted and observed, though shifted 
longitudinally by the same amount as the shocks. The two worst discrepancies 
involve the general underestimate of the temperature and the behavior of some 
of the other parameters, especially the density, following the reverse shock. 
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The former problem can be attributed to the single-fluid nature of the model 
(the T comparison is the one thing that can be significantly improved by 
recourse to the two-fluid formulation mentioned above), while the latter has 
its origins in the extended gap in the 1 AU data visible in Figure 4, left. 

In summary, then, the observed speed profile is accurately reproduced by 
the model, there i3 generally good agreement in the density and field strength 
profiles, and the temperature and pressure profiles show the same qualitative 
features. The relative positions of the forward fast shock and the stream 
interfaces are also accurately determined by the model. This is particularly 
significant because the shocks are handled by introducing a finite artificial 
viscosity to cope with rapid variations in plasma parameters. We conclude 
that the model is reliable over the distances considered and therefore feel 
justified in its further application in this paper. 

Finally, to illustrate how this stream relates to our preconceptions of 
CIR formation, Figure 7 presents a projection of the flow out to 4 AU. There 
the stream appears as the typical saw-tooth velocity structure, with forward 
and reverse shocks bounding and expanding compression region. The high speed 
material in the stream is being eroded by the reverse shock while the slow 
flow ahead of the CIR is being accelerated by the forward shock. The CIR 
continues to broaden with distance until eventually it begins to interact with 
neighboring CIR's (see also Burlaga, 1983; Pizzo, 1983). 
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5. Interaction of Adjacent Slow and Fast Corotating Streams 

In the Interval November 10-16, 1 977* IMP-8 at 1 AU observed one 
corotating stream flow Immediately followed by a faster corotating stream, as 
shown in a plot of 5 min averages in Figure 8. This compound stream is 
labeled "stream 1” in Figure 1, A stream interface was observed late on 
November 10, across which the density dropped while both the temperature and 
the bulk speed increased. As the speed continued to increase during the day 
following the interface, we may describe this structure as a slow corotating 
stream. There is some evidence for the formation of a forward shock on 
November 10 and possibly a reverse shock on November 11, which might be part 
of the same flow pattern. A relatively large and fast corotating stream was 
observed on November I 1 ! and 15, presumably originating in the trans-equatorial 
coronal hole shown in Figure 1. The magnetic field polarity of this stream 
differs from that of the preceding stream indicating that they had different 
sources. A diffuse interface between the two flows was observed on Novem- 
ber 14, though no forward or reverse shocks were seen by IMP-8. Because the 
corotating stream on November 14 and 15 was moving faster than that of the 
November 11 flow, one expects it to overtake, collide with and "sweep-up" the 
latter. The nature of that interaction is the principal topic of this Section. 

Figure 9 shows a plot of the connected hourly averages of B, V, N, T and 
P-j- (the magnetic plus proton pressure) versus relative heliocentric azimuth. 

The IMP-8 profiles on the left of Figure 9 were used as inputs to the MHD 
model, and the computed profiles at 2 AU are shown on the right of Figure 9. 

At 1 AU (Figure 9, left), both the magnetic field and the pressure show two 
maxima, one at each of the interfaces. Thus there are two compressions, one 
associated with the interaction region of the slow corotating flow, and the 
other associated with the interaction region of the fast corotating flow. The 
magnetic field profile is very similar to the total pressure profile, 
indicating the importance of the magnetic pressure in the dynamics. At 2 AU 
(Figure 9, right), in contrast to the two broadly spaced peaks in the pressure 
profile at 1 AU, there is a single peak. The same pattern is seen in the 
magnetic field profile: a single maximum in B is predicted at 2 AU whereas 

two widely-spaced maxima were observed at 1 AU. Thus, the model predicts a 
qualitative change in the total pressure and magnetic field profiles between 
1 AU and 2 AU. 



16 


The qualitative change just described occurs because the advection of the 
fast corotating stream relative to the slower flow brings their separate 
compression regions (also referred to in the literature as corotating pressure 
waves; Burlaga et al., 1983) closer together while at the same time each 
compression region wave is expanding. The amplitude of the stream at 2 AU 
(see Figure 9) is diminished beoause some of the kinetic energy of the 
high-speed material from the streams is continually being fed into the growth 
of their respective compression regions between 1 AU and 2 AU. The mechanism 
for this increase is kinematic steepening (e.g., see Hundhausen, 1972, Burlaga 
and Barouch, 1976; and Gosling, 1981). Thus the amplitude of the corotating 
pressure wave ahead of the fast stream grows relative to that ahead of the 
slower stream. At the same time, the expansion of the two compression regions 
tends to fill in and close the region of low pressure seen at 1 AU, so that by 
2 AU there is a single compression which is formed by the coalescence of two 
separate corotating pressure waves. 

The above scenario is qualitative and intuitive, of course, and while the 
validity of the model in the case of the simple corotating stream has been 
confirmed, it has not been unequivocally demonstrated under these 
circumstances. More quantitative results and additional insight can therefore 
be obtained by comparison of the predictions of the model with observations. 

In this case, we compare predictions of the model with data from Voyager 1 
that were obtained near 1.44 AU. The results are shown in Figure 10 with the 
predictions of the model shown by the solid curves and the Voyager 
observations shown by dots. 

Figure 10 shows that at 1.44 AU there were two broad humps in the 
pressure profile instead of the two narrower peaks seen at 1 AU. The model 
predicts the development of a reverse shock associated with the slow stream. 
This is the feature that was identified as a potential reverse shock on 
November 11 in Figure 8. A fully developed reverse shock was not seen by 
Voyager 1 but the observed general variation of the flow parameters near the 
predicted shock position did resemble the calculated pattern (see Figure 10). 
Thus, we may say at least that the interaction region ahead of the slow stream 
was probably expanding backward toward the sun at a speed comparable to the 
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magnetoaooustio speed, further broadening the interaction region, reducing the 
amplitude of the compression region, decelerating the stream, and filling the 
low pressure region observed at 1 AU. Of course, these features could not be 
derived from a kinematic model; dynamical effects are fundamental in the 
evolution of Interaction regions. 

The interaction region associated with the fast eorotating stream in 
Figure 10 shows a peak in near the interface, but the pressure wave is 
rather broad and flat-topped. Neither a corotating forward shock nor a 
reverse shock is observed or predicted, but one can see a region (marked (RS) 
in Figure 10) where a reverse shock might be forming. The forward part of the 
corotating pressure wave ahead of the fast stream is just approaching the 
reverse part of the pressure wave associated with the slow stream (which has 
the signature of a reverse shock in the theoretical curves in Figure 10). 

Thus, the low pressure region between the two pressure waves is narrow at 1.44 
AU and apparently being filled in rapidly. The agreement between predictions 
and observations in Figure 10 gives us further confidence in the accuracy of 
the model. In particular, the coalescence of the compression regions 
(pressure waves) predicted to occur between 1 AU and 2 AUmust be taken 
seriously, as the Voyager observations at 1.44 AU are covenlently located to 
reveal the process in mid-stride. We can now proceed to describe some of the 
other features predicted by the mode 1 . (Figure 9). 

Two reverse shocks, at 2 AU, on the sunward side of the interface, are 
predicted by the model. One of the reverse shocks is a corotating reverse 
shock associated with the fast stream while the other is the corotating 
reverse shock associated with the slow stream. The reverse shock associated 
with the slow stream moved through the interface associated with the fast 
stream before it reached 2 AU. This caused an additional compression of the 
corotating interaction region ahead of the fast stream, giving a further 
enhanced pressure at its interface. Of course, the low pressure region seen 
between the two pressure waves was also compressed by this reverse shock as it 
moved through the slow stream toward and past the interface of the fast 
stream. Thus, the observation of 2 reverse shocks may be a signature of the 
interaction between 2 corotating streams. The coalescence of the two pressure 
waves seen at 1 AU into a single pressure wave at 2 AU is seen as the result 
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of several processes; 1) the slow stream pressure wave expands and drops in 
amplitude; 2) the fast stream pressure wave is advected forwarded, but does 
not diminish, since it is 3till bfting driven by the fast material; 3) where 
they meet there is further compression. 

The model also suggests that other interesting interactions take place 
farther out, producing by 4 AU a single CIR with complex substructure. Figure 
11 shows the speed profiles and the magnetic profiles (which are similar to 
the pressure profiles, not shown) predicted at 1 AU, 2 AU, 3 AU, and 4 AU. 

The most significant qualitative change between 2 AU and 4 AU is the 
coalescence of the two reverse shocks near 3 AU. This occurs because the 
reverse shock associated with the slow stream was moving faster than the 
reverse shock associated with the fast stream. It also appears that a forward 
shock associated with the fast stream forms near 2 AU and overtakes the 
forward shock associated with the slow stream Just beyond 4 AU. 
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6 . Destruction of a Corotating Stream overtaking a Transient Flow 

Me shall now describe the overtaking of a slow transient flow by a faster 
corotating stream. Coalescence of the two interaction regions is again 
observed, and in addition the destruction of the corotating stream is 
observed. This event has been discussed by Burlaga et al, (1980), who 
described the evolution of the corotating stream and suggested that the 
disappearance of the corotating stream between 1 AU and 1,6 AU was a dynamical 
effect. We shall give further details of the observations and show that the 
MHD model supports the oonjecture that the corotating stream was destroyed by 
a dynamical interaction. 

The IMP-8 observations at 1 AU are shown in Figure 12 for the period 
November 25 to December 1, 1977. Three flow configurations can be seen: 1) a 

transient flow on November 25, consisting of a forward shock (FI), followed by 

-3 

a shell of very high density material (> 100 cm ), followed in turn by a 
region of strong magnetic fields moving relatively slowly; 2) a corotating 
stream with a well-defined stream interface (IF) on November 26, and 3) a 
shock (F3) on November 29. The shock F3 apparently did not interact with the 
other flows inside of 2 AU, so we shall first consider the interaction of the 
corotating stream with the transient flow ahead of it. 

Figure 13 shows the connected hour-averages of B, V, N, T, and P T measured 
at 1 AU together with the corresponding profiles at 2 AU computed from the MHD 
model. Ignoring for the moment the transient flow corresponding to F3» one 
sees two compression regions near 270° at 1 AU, but only one broad compression 
region at 2 AU. The two pressure waves at 1 AU correspond to the interaction 
regions of the transient flow and the corotating flow, respectively. 

The predictions of the model are compared with observations made at 1.6 AU 
by Voyager 1 in Figure 14. We see fairly good agreement between the 
theoretical speed profile and the observations, to within 20 km/sec over most 
of the time interval, and the basic qualitative features of N and B profiles 
are all reproduced by the model, with respectable quantitative agreement at 
higher values of N and B. The theoretical temperature and pressure profiles 
agree qualitatively with the observed profiles; there are some large 
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differences at times, which might be attributed to inhomogenelties of the 
flow, temporal effects, and to our neglect of electron measurements, but these 
are differences in detail. The position of the forward transient shock (FI) 
is accurately predicted. The model also prediots a corotating forward shock 
(F2) and a corotating reverse shock (RS), The corotating forward shook 
arrives slightly later than predicted, but this is again a detail which we 
shall not pursue in this paper. The model does accurately predict the time of 
arrival of the reverse shock. Detailed observations of the shocks themselves 
have been presented by Burlaga et al. (1980). 

Figure 13 shows a narrow corotating stream with a maximum speed of 
* 525 km/s associated with the interface at 1 AU, but at 2 AU the speed does 
not exceed 400 km/s; it is this change that we referred to above in speaking 
of the destruction of the stream. Since the 2 AU curves are theoretical 
profiles based on an HHD model of corotating streams, we have demonstrated 
that the destruction of a small corotating stream within 1 AU is theoretically 
possible for the input conditions given in Figures 12 and 13. (It should be 
emphasized that this case may not be typical, for the stream was relatively 
small and it collided with a region of unusually high density (n > 100/cm^) 
associated with the transient shock.) It is firmly established, of course, 
that streams are destroyed (in the sense that velocity differences are wiped 
out by dynamic interaction) over spatial scales of several AU (Hundhausen and 
Gosling, 1976; Collard et al., 1982). 

The two large pressure waves observed at 1 AU (see Figure 13 near 270°) 
are also seen at 1.6 AU (see the bottom of Figure 14). At 1.6 AU the 
corotating pressure wave was broader and had a higher amplitude than the 
transient pressure wave. It appears that these two waves were interacting and 
beginning to coalesce at 1.6 AU. The single broad pressure wave predicted at 
2 AU (see Figure 13) is a natural extrapolation of the result seen in Figure 
14. The corotating interaction region and the transient interaction region 
coalesce and become nearly indistinguishable by 2 AU, the result being a 
single merged interaction region. 

Let us now consider the evolution of these flows to 3 AU. Figure 15 
shows the speed and magnetic field strength profiles computed with the MHD 
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model, using as input conditions the IMP-8 data shown in Figure 13. One sees 
that the corotating forward shock F2 overtakes and Interacts with the 
transient shock FI at *4 AU, The corotating reverse shock RS and the 
transient shook F3 interact somewhere between 2 AU and 3 AU, At 3 AU one sees 
a single broad interaction region in the B profile (the corresponding pressure 
profiles are not shown but they are qualitatively the same as the B profiles). 
This ’’interaction region" actually represents the coalescence of 3 separate 
interaction regions—those associated with the transient shock FI, the 
corotating stream, and the transient shock F3, Thus, we have at 3 AU a single 
merged interaction region and this is not associated with a fast streaml 
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7. Summary 

We have discussed multispacecraft observations and an MHD model of the 
radial evolution of three flows between 1 AU and 1.74 AU using data from 
Voyager 1, IMP-7 and IMP-8. The first flow system (stream 3) was a simple, 
classic corotating stream. A forward-reverse shock pair was observed 
relatively close to 1 AU (v.z. 1,74 AU) , and its development was predicted by 
the model. There was a small error in the predicted position of the shocks, 
which might be due to the small latitudinal separation of the spacecraft 
and/or to a temporal variation in the stream a3 it rotated through the several 
degree angle corresponding to the azimuthal separation of the two spacecraft. 
The second flow system consisted of a fast corotating stream overtaking a 
slower corotating stream. Two compression regioins (pressure waves) were 
observed at 1 AU, and these were computed to coalesce into a single pressure 
wave at 2 AU. The process was observed at mid-point by Voyager 1, confirming 
the gross predictions of the model. The third flow system consisted of a 
corotating stream overtaking a transient shock and flow. Again we observed 
two pressure waves at 1 AU which coalesced to form a single pressure wave at 2 
AU. A third pressure wave associated with another transient shock at 1 AU was 
predicted to coalesce with the pressure waves just described, forming 't single 
pressure wave at 3 AU from the 3 distinct pressure waves observed at 1 AU. 

This event was particularly notable because the corotating stream was 
effectively destroyed by its interaction with the transient flow inside of 
v 1.6 AU. 

The flows just described were modeled very satisfactorily with a 
stationary, 2-D, single- fluid MHD model. Comparison of predicted and observed 
flow features in all three cases demonstrates the general quantitative as well 
as qualitative validity of such models, down even to such details as shock 
formation and inspire some confidence in their predictions of dynamical 
behavior outside the realm of observation. The applicability of the 2-D 
assumption was discussed by Pizzo (1983); the justification for the 1-fluid 
formulation was investigated by comparing the results of this model with those 
computed for a 2-fluid model based on separate polytropic relations for 
protons and electrons. This latter model is essentially identical to that 
described by Goldstein and Jokipii (1977)» the main difference being that ' 
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electron polytrope constant is chosen to be r 3 1.175. [See Scudder and 
Olbert (1979a »b) and Sittler and Scudder, 1980 for the theoretical and 
observational bases for using a polytrope description for solar wind 
electrons.] The differences in the V-T-N profiles for the two cases were 
negligible and the differences in the proton temperature profiles were small, 
Basically, this is because the contribution that the electrons make to the 
pressure is smaller than that due to the sum of the magnetic field pressure 
and the proton pressure, and hence the electron correction is negligible. The 
assumption of stationarity is valid for corotating structures like stream 1 
and stream 3, but it is not obvious that it is appropriate for cases like 
stream 2, which involves a transient ''s well as a corotating flow. In fact, 
the model could not be applied if the spacecraft were not radially aligned. 

If the azimuthal separation were suitably large the transient flow might not 
even be seen by both spacecraft! However, we have analyzed a situation in 
which the spacecraft were nearly radially aligned, and we considered the 
evolution over relatively small distances. Under these circumstances , details 
of the azimuthal variations of the transient flow are not important to first 
approximation, since the effects of radial gradients are included in the 
stationary model. Thus, the good agreement between the observations and the 
predictions of the model is not Just coincidence, but rather stems from the 
fact that we observed mainly radial variations rather than a mixture of radial 
and azimuthal variations. 


8. Discussion 

The above results provide a graphic and compelling demonstration of how 
the magnetic field and plasma patterns can change as solar wind flows interact 
between 1 AU and a few AU. On the one hand, shock formation in a simple, 
isolated corotating stream does appear to take place just as qualitative and 
quantitative theory would suggest. On the other hand, in the more common 
scenario of a solar wind with neighboring, interacting corotating and 
transient flows, it is seen that structures in the more distant solar wind, 
say at 4 or 5 AU, will be in fact a complex amalgam of dynamically merged 
flows, even though at that distance they may bear at least superficial 
resemblance to text-book CIR T s resulting from large, isolated streams at 1 AU. 



The dynamical process by which adjaoent flows of varying spatial and 
temporal scales merge together at larger heliocentric distances has been 
described in the context of filtering (Gosling et al., 1976; see also the 
review by Holzer, 1979). That i3, one of the most fundamental dynamical 
properties of the solar wind is the tendency for small scale flow structures 
to evolve more rapidly than large scale structures, such that at greater 
heliocentric distances only the largest scale structures remain, the smaller 
one3 having been washed out (through dynamical exchange of momentum and energy 
with their surroundings) and subsequently 3wept up by the surviving larger 
one3. It is in this sense, then, that the 3olar wind is said to act as a 
low-pass filter. A fairly straightforward manifestation of this prooess is 
evident in our stream It 3, the classic isolated flow case. Referring back to 
Figures 4 and 5, we call attention to the high frequency (short wavelength) 
fluctuations discerned in the 1 AU IMP velocity profile. By 1.74 AU (Figure 
5), both the calculation and the real data show a marked loss of detail, a 
distinct snoothing of the profile, The effeot is particularly noticeable in 
the trailing portions of the stream. Gosling et al. (1976), remarking upon 
the same phenomenon in their 1 AU to 4 AU mapping of IMP-Pioneer streams, 
ascribed the disappearance of the fluctuations in speed to rapid (relative to 
that of the overall waveform) steepening commensurate with their small scale, 
leading to shock formation andd degradation via dissipation. While shock 
formation and dissipation may have been a significant ingredient of the 
proaes3 within the context of the idealization of their model (an r-t 
hydrodynamic formulation known to overestimate the intensity of such 
interactions by virtue of its neglect of magnetic and nonraa'ial flow effects), 
it is really only incidental to the dynamical evolution. That is, the 
essential process is that small scale structures evolve more rapidly and are 
eventually consumed by larger scale features. If the amplitudes of the small 
scale features are sufficiently large, shocks may form and indeed hasten the 
demise of the fluctuations; however that is not at all essential to the 
process, which can proceed to completion (albeit perhaps more slowly) by 
following a ’'subsonic" shock- free evolution (Hundhausen and Pizzo, unpublished 
manuscript; see also the discussion of "irregular fluctuations" in Burlaga, 
1975). Finally, it should also be emphasized that it makes no fundamental 
difference whether the small scale structures are corotating or truly 
temporally varying; the evolution may proceed at different rates, but it will 
occur all the same, as demonstrated by streams 1 and 2 in this study. 
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For the smallest scale, smallest amplitude fluctuations, the dynomioa/ 
evolution pretty much leads to oblivion, in the sense that all traces of their 
intrinsic velocity, thermodynamic, and magnetic structure are ultimately 
washed out beyond recognition by the larger scale, higher amplitude flows in 
which they beoome embedded. For the larger fluctuations, however, it is to be 
expected that their intrinsic structure might have an important Impact upon 
the process and that some residue of the original variations might even be 
disoernible at laarger heliocentric distances. This aspect of the 
evolutionary dynamics was also addressed in Gosling et al. (1976). In their 
examples, the individual velocity variations associated with several umall 
scale streams lying at the leading edge of a muoh larger stream were observed 
to be smoothed out and virtually obliterated by 4 Al); but their associated 
compression fronts (which were quite distinct at 1 AU) had been merged 
together into a dynamical evolution. The mechanism for this merging, of 
aourse, is simply the kinematic advection of the larger scale, higher 
amplitude streams into the smaller ones together with the interaction of the 
individual pressure waves. Unfortunately, the Gosling et al. (1976) study was 
somewhat handicapped, in the sense that the start (1 AU) and end (4 AU) points 
of the process were all that could be observed (and then only the velocity 
profiles were available at 4 AU); the authors had to rely upon intuition and a 
nunerical model to reconstruct what had happened in the intervening distances. 
The contribution of our study to this matter is that the IMP and Voyager data 
utilised herein, which are rather complete, captures the process in action and 
it shows that even quantitative details of the evolutionary process can be 
tracked with the appropriate model. 

Recently, the importance of the cummulative effects of the dynamical 
interactions among various corotating and transient flows in the inner solar 
system upon the properties of the distant solar wind has received special 
emphasis (Burlaga, 1983). Indeed, the specific sweeping and merging (or 
coalescence) aspects of stream evolution have been accorded the specific 
designation "entrainment" (Burlage et al., 1983a). The entrainment of 
compression regions (such as occurs in our streams //I and it2) has been shown 
to imply the transfer of magnetic energy from smaller to larger scales 
(Burlaga 'and Goldstein, 1983). It has been suggested that the magnetic field 
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structure of various kinds of flow systems (e.g., those consisting of 
predominantly corotating structures versus predominantly transient structures) 
at large heliooentrio distances might differ somewhat with significant impaot 
upon cosmic ray transport (Goldstein at al., 1984} Burlaga et al., 1983b). 

Finally, the dynamical filtering of small scale, stream structures may 
also be taking plaoe well inside 1 All, though perhaps with some modifications. 
This is, near the sun, stationary corotating fluctuations could be expected to 
evolve only slowly, as the rotational mechanism driving the dynamios produces 
a very oblique interaction between adjaoent stream tubes. Non-stationary 
fluctuations close to the sun, on the other hand, would in general possess 
more nearly radial gradients, and we might anticipate that they would be 
subjeot to stronger interactions. It is interesting to speculate that this 
tendency for the temporal fluctuations to be more strongly filtered near the 
sun which might have some bearing upon the afore-mentioned difficulties 
encountered by the models in treating the near sun evolution of cororating 
streams (see also Pizzo, 1983, 1984). Thi3 is, considerable low-amplitude 
fluctuations associated with the ever-present global coronal evolution should 
be present in the solar wind in the inner solar system, (Here we expressly 
exclude consideration of impulsive, large scale, transient systems, such as 
associated with flares.) Hence the evolution cannot strictly follow that 
expected for a single, absolutely corotating stream. Near and beyond 1 All 
however, nearly all these short-period temporal variations would have been 
washed out, with only the largest amplitude, longest wavelength structures 
remaining. In this view, then, the relative success of corotating stream 
models at larger heliocentric distances cannot be attributed to the geometric 
diminution of nonradial flow effects alone. Rather, flow interactions and 
filtering may play a significant role, in that they act in such a way as to 
assure that the greater the heliocentric distance, the closer the remnant flow 
systems oome to the steady idealization presumed in the models. 
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(Top) Speed profile measured by Imps 7 and 8 showing the three 
flows discussed in the paper; the polarity of the inter- 
planetary magnetic field is also shown. (Bottom) K-coronameter 
polarization brightness contours, showing the location of 
coronal holes. 

Trajectories of Imps 7, 8 and Voyagers 1, 2, in inertial 
heliographic coordinates. The angle X is the heliographic 
latitude of the spacecraft. 

Five minute averages of the magnetic field and plasma data for 
stream 3, showing a stream interface at the leading edge of a 
corotating stream. 

(Left) Magnetic field strength, plasma parameters, and total 
pressure for stream 3 observed at 1 AU. (Right) The 
corresponding parameters calculated for 2 AU, showing a forward 
shock (FS) and a reverse shock (RS). 

Magnetic field strength, plasma parameters and total pressure 
versus Carrington longitude. Predicted values from the MHD 
model with IMP data at 1 AU as input are shown by the solid 
curves. Observed values are shown by the dots. 

High temporal resolution data of stream 3 from Voyager 1, 
confirming the existence of a forward-reverse shock pair at 
1.74 AU. 

Theoretical curves for the speed and field strength at 2, 3 and 
4 AU, together with the profiles measured at 1 AU. 

Five minute averages of the magnetic field and plasma data for 
stream 1, showing a forward shock (FS) , two stream interfaces 
(IF), and a possible reverse shock (RS?). 
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Figure 9 (Left) Magnetio field strength , plasma parameters and the total 

pressure (magnetic plus thermal) for stream 1 observed at 1 AU. 
Two compression regions are indicated by the shaded areas in 
P.p one for the slow stream and one for the fast stream. 

(Right) the c.or responding parameters calculated for 2 AU. A 
single broad compression is visible at 2 AU, implying the 
coalescence of the two shown on the left. 


Figure 10 Magnetic field strength, plasma parameters and total pressure 

versus Carrington longitude. Predicted values from the MHD 
model with IMP data at 1 AU as input are shown by the solid 
curves. Observed values are shown by the dots. 

Figure 1 1 Theoretical curves for the speed and field strength at 2, 3, 

and 4 AU, together with the profiles measured at 1 AU. Note 

the coalescence at the two peaks in B between 1 AU and 2 AU, 

the coalescence of the 2 forward shocks near 4 AU, and the 

coalescence of the 2 reverse shocks near 3 AU, 


Figure 12 Five minute averages of the magnetic field and plasma data for 

stream 2, showing a forward shock (FI), an interface (IF), a 
possible site for formation of a reverse shock (RS), and a 
possible forward shock (FS). 

Figure 13 (Left) Magnetic field strength, plasma parameters, and total 

pressure for stream 2 observed at 1 AU. (Right) The 
corresponding parameters computed for 2 AU. Two pressure waves 
observed at 1 AU coalesce to form a single pressure wave at 
2 AU. 


Figure 14 Magnetic field strength, plasma parameters and total pressure 

versus Carrington longitude. Predicted values from the MHD 
model with IMP data at 1 AU as input are shown by the solid 
curves. Observed values are shown by the dots. 
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Figure 15 Theoretical curves for the speed and field strength at 2 AU and 
3 AU, together with the profiles measured at 1 AU. Note the 
coalescence of the shocks and of the pressure waves. 
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